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PROCESS AND APPARATUS FOR DETECTION OF THE CARDIAC CYCLE 
CROSS -REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims a benefit of a priority under 35 USC 1 19(a)-(d) to 

FR 02 15443 filed December 6, 2002, the entire contents of which are hereby 
incorporated by reference. 

BACKGROUND OF THE INVENTION 

[0002] This invention relates to a process and an apparatus for detection of a 

cardiac movement or cycle of a patient. More particularly, the invention relates to the 
detection of the cardiac cycle starting from an angiogram of coronary vessels. 
Identification of cardiac movement from images of coronary arteries acquired using a 
radiography apparatus during selective injection of a contrast agent within these coronary 
arteries, can be useful in a number of algorithms: 

[0003] Calculation of a myocardial perfusion index based on analysis of the 

variation of density in a series of images throughout different cardiac cycles. 

[0004] Production of a fixed sequence, in other words a sequence with no or very 

few displacements of the coronary arteries between the different images making up the 
sequence. 

[0005] Synchronization of two multi-modal sequences for the same patient, for 

example, one sequence acquired in fluoroscopy and another sequence acquired in 
conventional recording mode, the images of the two sequences possibly being combined. 

[0006] Production of a three-dimensional image reconstruction by applying 

known reconstruction algorithms from a series of images chosen to illustrate the same 
instant in the cardiac cycle. 
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[0007] There are difficulties with these algorithms. A first difficulty is that the 

electrocardiogram (ECG) signal is necessary to know the instant in the cardiac cycle 
associated with each image acquired in the sequences. But this signal is not always 
directly available with the images. Furthermore, the ECG signal is an electrical signal 
that describes the stimulation of the heart muscle but not the effective mechanical 
movement of this muscle. However in the known application of algorithms, it is 
necessary to know the position of injected coronary arteries, in other words to know the 
mechanical behaviour of the heart 

[0008] A second difficulty is that the patient might breath during the acquisition 

or that the images from two different cardiac cycles may be slightly different. 

BRIEF DESCRIPTION OF THE INVENTION 

[0009] An embodiment of the invention provides a process for the detection of 

the cardiac movement capable of solving the above-mentioned difficulties. An 
embodiment of the invention, is a process for a cardiac movement detection designed to 
be implemented within a radiography apparatus of the type comprising means for 
providing a source of radiation, means for recording placed facing the source, and a 
means for support placed between the source and means for recording on which there is a 
patient for whom an image of a region of the heart is to be made. The process may be 
automatic and implemented by a suitable computer program. 



[0010] An embodiment of the process comprises: 

[0011] a) acquisition of a series of successive images I n of a region of a 
heart; 

[0012] b) analysis of at least some of the images thus acquired to identify a 



heart movement; and 
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[0013] c) determination of the cardiac cycle starting from this movement. 

[0014] Thus, an embodiment of the process can be used to determine the effective 

mechanical movement of the heart during its cycle starting from a series of images 
acquired by a radiography apparatus without the need for the electrocardiogram signal so 
that the different images can be synchronized with the electrocardiogram. Furthermore, 
during the calculation of integral displacements between successive images in the series 
of images, an embodiment of the process can evaluate minor differences that may arise 
between two cardiac cycles contained in the series of images. 

[0015] An embodiment of the invention is a radiography apparatus of the type 

comprising means for providing a source of radiation, means for recording facing the 
source and means for support placed between the source and the means for recording on 
which there is a patient for whom a region of the heart is to be imaged, the radiography 
apparatus comprising means for implementing the process. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0016] Other characteristics and advantages will be better understood from the 

following description of an embodiment of the invention and the drawings, in which: 

[0017] Figure 1 is a diagrammatic view of a radiography apparatus capable of 

implementing the process; 

[0018] Figure 2 is a principle diagram of the process; 

[0019] Figure 3 shows two successive images acquired during the process; 

[0020] Figure 4 is the transformation of the images in figure 3 into segmented 

images; 
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[0021] Figure 5 is a diagram representing the attenuation coefficients calculated 

using the process for the images in Figure 4; and 

[0022] Figure 6 is a diagram illustrating the cardiac cycle obtained using the 

process. 

DETAILED DESCRIPTION OF THE INVENTION 

[0023] With reference to Figure 1, a radiography apparatus 1 comprises means for 

recording images 2 and means for providing a source of radiation 3, such as an X-ray 
source. The means for recording images 2 may comprise radiographic plates of a plane 
sensor or a luminance amplifier associated with a camera. The radiation source 3 and the 
means for recording images 2 are fixed to each end of a carrier arm 7 acting as a 
positioner, in this case in the form of a semi-circle. The semi-circular arm 7 is connected 
by sliding on a second arm 8. The second arm 8 is connected in rotation about an 
installed base 9 of the radiography apparatus 1. The base 9 is free to rotate 12 with 
respect to the ground. 

[0024] Arm 8 is essentially capable of making rotation movements 6 about its 

own axis. The semi-circular arm 7 is capable of sliding with respect to the arm 8, such 
that the semi-circular arm 7 makes a rotation movement 5 with respect to the center of the 
semi-circle forming the arm 7. 

[0025] During use, the body of the patient 100 is positioned between the radiation 

source 3 and the means for recording images 2, on a means for support (not shown) such 
that the heart 4 of the patient 100 is within the field 10 of the apparatus. 

[0026] With reference to Figure 2, the patient 100 is lying down on the means for 

support of the radiography apparatus 1 while the radiation source 3 and the means for 
recording images 2 make a rotation 5 around the patient. In this case, the cranio-caudal 
axis of the patient 100 is approximately parallel to the axis of rotation 5 of the movement. 
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Anatomically, during the cardiac cycle, two planes of the heart perpendicular to its 
cranio-caudal axis make a relative translation movement along this cranio-caudal axis 
associated with a rotation movement about this cranio-caudal axis. Consequently, with 
the patient 100 being positioned within the radiographic apparatus 1, each plane P t of the 
heart perpendicular to the cranio-caudal axis of the heart is projected along a line Lt onto 
the means for recording images when this means acquires an image S t . 

[0027] At time t+5t, when the semi-circular arm 7 supporting the radiation source 

3 and the means for recording images 2 has made a rotation movement 5 about its axis, 
the previously described plane P t has made a translation about the cranio-caudal axis of 
the heart and a rotation about this same axis and is in the position P t +8t that is then 
projected onto means for recording images 2 along a line Lt+5 t that is shown on image 
S t +st acquired at this moment. The P t and P t +8t planes intersect the same types of cardiac 
tissues in the same proportions. Since tissues attenuate radiation passing through them by 
variable amounts, the sum of these attenuations along line Lt is approximately equal to 
the sum of the attenuations along line Lt+8t. Thus, to determine the displacement of plane 
P, all that is necessary is to calculate the displacement of its projection on the means for 
recording images 2, in other words the displacement of line L from Lt to Lt+§ t . 

[0028] Consider Roo, the intensity of the ray emitted by a radiation source 3. The 

image I obtained after injection of a contrast product in the coronary vessels of the heart 
may be modeled by: 

[0029] I(M) = Rooe c ^ {x)dx ~i c ^ x)dx 

[0030] where C(M) is the path between the radiation source 3 modeled in Figure 

2 by point O and a point M on the image, the point M having coordinates (i J) in pixels on 
the image; ji is the local attenuation coefficient along the path C that depends on the 
nature of tissues passed through and the wavelength of the X-rays used; V represents all 
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points on the image belonging to the projected vessels through which the radiation 
passes; and F represents all points belonging to other tissues projected on the image. 

[0031] If the image were not injected, the contributions of points that are not in 

vessels would be obtained, in other words an image Io that would be modelled as: 

[0032] h(M) = Rooe-l c(F)fi{x)dx 

[0033] Consequently, the principle of the process is based on the fact that 

J Cin M(x)dx « s conserve( j Tkjg en tity is calculated from different images in a series of images 
and the vertical displacement of the heart between two successive images in the series of 
images is then determined. 

[0034] In one embodiment, the images acquired by the means for acquisition, i.e., 

the radiographic plates 2, are images of the type as shown in Figure 3. In Figure 3 an 
image IO n and the next image IO n +i in the series, on which the later images are visible 
and dark due to the injection of contrast products into the coronary vessels. The first step 
in the process, knowing the maximum dimensions of the arteries to be extracted, in other 
words a diameter of the coronary arteries, is to apply a closing mask to the initial images 
IO n as illustrated in Figure 3. This closing mask may bring out structuring elements with 
a size equivalent to the dimensions of the arteries to be extracted and determined in 
advance. Application of this type of mask involves calculating a subtracted image 
obtained by subtraction of the closing mask from the corresponding previously calculated 
initial image. However, this operation does not eliminate all parasite structures present in 
the image. In order to achieve this, the process applies adaptive thresholding to the 
subtracted image by calculating a threshold so as to keep about 1 5% of the pixels. This 
assumption is based on the fact that the coronary vessels that contain a contrast product 
when the images are acquired give the greatest contrast. The result is a thresholded 
contrast image like that illustrated in Figure 4 in which the image IS n is the thresholded 
image of the original image IO n , and similarly, the thresholded image IS n +i is the 
thresholded image of the original image IO„+i. 
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[0035] In a second step, the process determines the linear integrals f ISu of the 

different thresholded images IS n , pixel line i by pixel line i (in this case these pixel lines 
are perpendicular to the cranio-caudal axis of the heart), according to the formula: 

No. of columns 

>*<•> &»• 

j=0 

[0036] The result is that i values of linear integrals f ISn illustrated in Figure 5 are 

obtained for each image IS n , where IA n is the set of linear integrals associated with the 
image IS n and IA n +i is the set of linear integrals associated with the image IS n +i. 

[0037] In the next step, the process calculates an integral displacement on a 

sequence of images starting from all previously calculated linear integrals IA n and 
associated with each original image IO n . The calculation of the integral displacement on a 
sequence of images determines the effective mechanical movement of the heart, or at 
least determines the variation of the translation of the different planes perpendicular to 
the cranio-caudal axis of the heart along this axis, and for which the projection of each of 
these planes on the images IO n is a pixel line i. In order to achieve this, the previous 
description indicated that the sum of the attenuation coefficients of projection points from 
the plane P t perpendicular to the cranio-caudal axis is conserved, in other words the 
values f ISn of the linear integrals are conserved. The process calculates the displacement 

between two successive images IO n and IO n +i along this translation. The process then 
determines a displacement kn, n +i starting from linear integrals IA n , IA n +i associated with 
each of the successive two images considered, k n , n +i that minimizes a cost function F n , n +i 
defined as follows: 

No. of lines . 

[0038] F n , n+1 (k n , n+ ,)= £ |/ /s>i (i)-/ /Vi (i-k nn+l) | 

i=0 

[0039] The value k n , n +i thus found is equal to the effective translation 

displacement of the heart along its cranio-caudal axis between two successive images 
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considered. The set of values kn, n +i found for the entire image sequence shows the cardiac 
cycle as illustrated in Figure 6. 

[0040] According to an embodiment, in an additional step, the process assumes 

that the sequence of events is approximately the same in each cardiac cycle, and 
determines the synchronous images of an arbitrary origin image taken from the first 
previously found cycle. However, in another variant embodiment for which all 
synchronous images to be found should enable a three-dimensional image reconstruction, 
images for which the contrast product is actually injected into the coronary vessels taken 
at the beginning of rotation of the semi-circular arm 7 of the radiography apparatus 1 
should be taken into account. Different tests have shown that approximately one hundred 
well-injected images can be used out of a sequence containing a total of about two 
hundred images, and at least three synchronous images can be found among the one 
hundred well-injected images that can be used for three-dimensional reconstruction. 



[0041] In one variant embodiment, the process calculates displacement of the 

heart due to the breathing movement of the patient during acquisition of the sequence of 
images used as input data to the process. In this step, the process is based on all 
previously chosen synchronous images. In a manner very similar to what was presented 
above, the process calculates the heart displacement due to the breathing movement. In 
order to illustrate this movement, it is assumed that the process uses three synchronous 
images Ii, I2, 13 and associated thresholded images ISi, IS2, IS3. Starting from these three 
images, a cost function F123 is defined that takes account of the fact that translation 
between images Ii and I3 is equal to the sum of the translation between images Ii and I2 
and the translation between images I2 and I3, in other words k^ = kn + k23. 

[0042] The expression of this new cost function F123 is: 

No. of lines . No. of lines . . No. of l ines 

FI 1 23(k 1 2,k 23 )= £ l/ JS) (iK/s/^,2> >+ £ l/^GHis^'-^l + Z 

i=0 i=0 i=0 
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[0043] Thus, the displacements ki 2 between images Ii and I2 and IC23 between 

images I 2 and I3 indicating displacements due to the breathing movements between image 
Ii and image I2 and image I2 and image I3 respectively, where ki2 and K23 minimize F123. 
Thus, before making a three-dimensional image reconstruction starting from synchronous 
images, the radiography apparatus 1 corrects these images starting from displacements 
due to breathing previously calculated by known image correction algorithms by 
translation, and then by known projection matrix correction algorithms. 

[0044] Depending on the variant embodiments, the process may comprise at least 

one of the following characteristics: 

[0045] (a) acquisition of a series of successive images I n of the region of the 

heart; 

[0046] (b) determination of a cranio-caudal axis of the heart; 

[0047] (c) for each image I n , calculation of the series of images of a set of 

attenuation coefficients of points on the image representing vessels in the region of the 
heart along lines perpendicular to the cranio-caudal axis of the heart; 

[0048] (d) calculation of the integral overrun k n , n +i between two successive 

images I n and I n +i of the series of images starting from the set of attenuation coefficients 
calculated for each of the two successive images; and 

[0049] (e) determination of the cardiac cycle starting from all previously 

calculated integral displacements. 

[0050] The determination of a cranio-caudal axis of the heart may comprises: 
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[0051] (bl) for each image I n in the series of images, calculation of an associated 

thresholded image IS n , only keeping vessels in the heart region; and 

[0052] (b2) determination of all attenuation coefficients for points on the image 

starting from thresholded images along lines perpendicular to the axis of the heart. 

[0053] The calculation of the thresholded image IS n may comprises: 

[0054] (bli) determination of at least one dimension, particularly a diameter, of 

vessels in the heart region to be kept; 

[0055] (blii) calculation of a closing image starting from the maximum 

dimension of vessels in the heart region to be kept; 

[0056] (bliii) calculation of an intermediate image by subtracting the closing 

image from the initial image; and 

[0057] (bliv) calculation of the thresholded image by application of an 

appropriate thresholding on the intermediate image. 

[0058] Thresholding is applied such that it keeps only about 15% of the pixels in 

the intermediate image; 

[0059] the set of attenuation coefficients of points on image I n along line i is 

modeled by a linear integral /in(i) of these attenuation coefficients along this line; 

[0060] the linear integral is expressed by a formula for line i, /in(i) = 

No. of columns 

Yin^H where l n (ij) = Rooe C( ">" M *" s c ^" (x)dx and Io„(ij) = Rooe- \ cin Kx)dx, 

JmO " 

where Roo is the initial intensity of radiation emitted by the radiation source; C(M) is the 
path between the radiation source and the point M on the image I n with coordinates (i j) 
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in pixels on the image; is the local attenuation coefficient along path C that depends on 
the nature of the tissues crossed and the wavelength of the radiation used; V represents all 
points on image I n belonging to the projected vessels through which the radiation pass; 
and F represents all points belonging to other tissues projected onto image I n . 

[0061] The integral displacement kn, n +i between two successive images I n and I n +i 

is calculated starting from all linear integrals associated with each successive image; and 

[0062] the integral displacement is the value of k n , n +i that minimizes an 

F n , n +i(kn,n+i) = £ I f h (0"/,^ (i-k w>n+1 ) I type cost function. 

i 

[0063] The process may comprise an additional step (f) to choose a subset of 

synchronous images in the heart cycle from the series of images, starting from the 
previously determined cardiac cycle. 

[0064] The process may also comprises a step (g) to determine an integral 

displacement due to breathing of the patient (100) between synchronous images, the 
determination being done in the same way as in the step (c). 

[0065] One skilled in the art may make or propose a modification in 

structure/way and/or function and/or result for the disclosed embodiments and 
equivalents thereof without departing from the scope and extant of the invention. 
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